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Abstract 

We determine the cc component of the photon wave function on the basis of 

(i) the data on the transitions e+e" -> J/V(3096), -0(3686), ^(4040), ^(4415), 

(ii) partial widths of the two-photon decays r] C Q (2979), XcO (3415), Xc2 (3556) — > 77, and 
(hi) wave functions of the charmonium states obtained by solving the Bethe-Salpeter 
equation for the cc system. Using the obtained cc component of the photon wave func- 
tion we calculate the 77 decay partial widths for radial excitation 2S state, r/ c o(3594) — ► 
77, and 2P states XcO (3849), Xc2 (3950) -► 77. 

1 Introduction 

There is a number of processes that can be sensibly treated by introducing the cc component 
of the photon wave function. First, it is the production of charmonium in the two-photon 
transitions such as 7*7* — > cc-mesons, production of -^-mesons in the e + e~-annihilation and 
production of charmonia in photon-nucleon collisions, 7*p — > cc-meson + X. In the present 
paper, we determine the cc-component of photon wave function, or the transition vertex 
7 — > cc, following the method developed in [1,2], where quark-antiquark components were 
found for the transitions 7 — > uu, dd, ss. 

The method of introducing quark-antiquark photon wave function may be clearly illus- 
trated by considering the cc-meson — > 77 decay. Dealing with the time-ordered processes, 
that is necessary in the dispersion relation approach or light-cone variable technique, the 
cc-meson — > 77 decay is a two-step process: first, the emission of photon by quark (Fig. la) 
or antiquark (Fig. lb) and, second, subsequent annihilation cc — > 7. 

In Refs. [1, 2], the triangle diagrams of Figs. la,b were treated in terms of double 
dispersion relation representation. The double spectral integral cuttings of the diagram in 
Fig. la are shown on Fig. lc. In the diagram of Fig. lc, on the left from the first cutting, 
there is the transition vertex of charmonium, G c h arm0 nium(s), decaying into cc pair, where s 
is the quark invariant energy squared. In light-cone variables 

x(l — x) 
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where m c is the mass of c-quark and (x, kj_) are the characteristics of one of the quarks (frac- 
tion of the momentum along the z-axis and transverse component). In the dispersion integral, 
the left-hand cutting leads to the factor G charmonium / (s - M 2 charmon ium ) , where M charmonium 
is the charmonium mass, this factor being the wave function of the initial cc-meson: 

G charmonium\S) t / \ /r>\ 

™ charmonium\S ) ■ 



s _ A/f 2 

charmonium 



Likewise, the right-hand cut in Fig. lc, by describing the transition cc — > 7, gives us the 
factor: 

-e c , (3) 

where s' is the invariant energy square of quarks in the final state and e c is the c-quark 
charge. 

When we deal with the transition cc — > 7, the interaction of quarks should be specially 
taken into consideration. The quarks may interact both in the initial (Fig. 2a) and final (Fig. 
2b) states. In fact, the interaction of quarks in the initial state has been accounted for by Eq. 
(2), because the vertex function G c h armon i um (or wave function ^charmonium) is the solution 
of Bethe-Salpeter equation — diagrammatically, this equation is shown in Fig. 3a. As to 
quark interaction in the final state, it should be particularly taken into account in addition 
to the point-like interaction (3). The diagram shown in Fig. 3b stands for the description 
of quark interaction in the transition cc — > 7, and we approximate it with the sum of the 
■^-mesons pole terms, see Fig. 3c. Accordingly, the factor related to the right-hand cut of 
Fig. lc is written as follows: 

G ^f S>) e c , (4) 

where the vertex function G 7 _^ cg (s') at s' ~ 4m 2 c is a superposition of vertices of the S'-wave 
■?/>-mesons (see Fig. 3c): 

G^ c5 (s) ~ ^C n G^ nS ){s) , s ~ Am 2 c ■ ( 5 ) 

n 

Here n is the radial quantum number of ip- meson and C n 's are numerical coefficients — their 
definition is the task of this paper. At large s, the vertex cc — > 7 is a point-like one: 

G^cc(s) ~ 1 at s > s • (6) 

The parameter sq can be determined using the data on e + e~-annihilation into hadrons: it is 
defined by the energy range, where the ratio R(s) = cr(e + e~ — > hadrons)/(x(e + e~ — > fi + fi~) 
reaches a constant-behaviour regime above the threshould of the charm production, R(s) ~ 
10/3. The data support the value s ~ (10 - 15) GeV 2 [3]. 

Therefore, to describe the transition cc — > 7 we may introduce a characteristics, which, 
similarly to (2), may be called the charmed quark component of photon wave function: 

%^£l = * T(?2) _-( S ) , (7) 

ft q 
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here q is the photon four-momentum. Let us emphasize that such a wave function is deter- 
mined at s > Am 2 c . 

There exists a reaction which is directly related to the photon wave function. This is the 
transition e + e~ — > ip(nS), see Fig. 4: here the loop diagram is defined by the convolution 
of meson wave function and the vertex 7 — > cc: ^/ n s <E> G 7 ^ cg . 

Dealing with the cc interaction in the transition 7 — > cc, we take into consideration the 
S-wave -?/>-mesons only, while the contribution of mesons dominated by the D-wave such as 
■^(3770) and ■0(4160), is neglected. The error coming from such a neglect can be evaluated 
by comparing the -0-meson production cross sections for S- and D-wave states in the e + e~- 
annihilation: it is of the order of 10%. The -D-wave admixture into the low-lying 1 -mesons 
is also small: it is of the order of 1% for J/ip(lS) and ip(2S) [4]. 

We are performing calculations of two- and three-point loop diagram in the spectral 
integration technique. All the equations used, up to trivial substitutions of quark masses 
and charges, were obtained in [5]. Because of that, in this paper we present final expressions 
only accompanied by necessary comments and references. 

The main difference between calculations for 7 — > cc and those for 7 — > uu, dd, ss carried 
out in [1, 2] consists in our regard towards quark wave functions. In [1, 2], we used phe- 
nomenological quark wave functions for 7T°, 77, 77', while for light vector mesons (p,u,(f>) we 
assumed the quark wave functions to be similar to analogous pseudoscalar-meson wave func- 
tions, with whom they form the lowest 36-plet in terms of SU(6)-symmetry. But here, when 
reconstructing the wave function for the transition 7 — > cc, we have applied the charmonium 
wave functions found out from the solution of Bethe-Salpeter equation [4] . 

It was a long history of calculation of charmonium states, and now there is a rich collection 
of results obtained in the framework of the nonrelativistic approaches [6, 7, 8, 9] as well 
as in different variants of the Bethe-Salpeter equation [10, 11, 12, 13], see also references 
therein. However, one should bear in mind that in calculations related to quark-ant iquark 
systems, one focuses as a rule on the description of levels (or masses) of the system. Had 
the potential or its relativistic analogue been known, the Bethe-Salpeter equation would 
undoubtedly provide us with both levels and wave functions. But the problem is that in 
fact our knowledge about the quark-quark interaction in the soft region is rather poor. 
Therefore, in the reconstruction of quark-antiquark levels one can obtain a variety of wave 
functions. Within the spectral integration technique used here, we see that unambigous 
determination of interaction is possible in the simultaneous description of both levels and 
wave functions, see discussion in Section 2.6.3 of [5]. Because of that, by describing the cc 
system, we have used in [4] as an input both the known levels and known values of radiative 
transitions (cc) in — > 7 + (cc) out . These radiative transitions are rather sensitive to wave 
functions, thence a comparatively good description of radiation transitions in [4] allows us 
to believe that the wave functions of lowest cc states are determined reliably too. 

The paper is organized as follows. In Section 2 we present briefly the information about 
the spectral-integration Bethe-Salpeter equation, in the framework of which the description 
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of the cc systems was performed in [4]. In Section 2 of the present article we also give the 
charmonium wave functions obtained at the determination of G^ cS . In Section 3 we write 
down the formulae for the transition amplitudes e + e~ — > ip and r) c0 , Xco, Xc2 77, which 
are used in the fit, and present the results for the photon wave function. Brief summary is 
given in Conclusion. 



2 Charmonium wave functions 

The spectral integral equation for quark- ant iquark wave function, which can be convention- 
ally called as the Bethe-Salpeter equation, is written for a system with the total momentum 
J, angular momentum L, quark-ant iquark spin S and radial number n. For the cc system 
it reads: 

( S ~~ ^charmonium) *(n) /li-fij (^-0 = (8) 

= J ^d® 2 {P'-XX)V{s,s\{k ± k> ± )){k[ + m . 

Here the quarks are mass-on-shell, k\ = kf —k\ — k 2 = ml, and the phase space factor in 
the intermediate state is determined as follows: 



d<S> 2 (P'; k[, hQ = f*' 1 ,, , ^ft- {2ix) A 5 {A \P' -k[- k' 2 ) . 
2V ' x ' 2J 2( y 2n) 3 2k' w ( y 2n) 3 2k' 20 K ; v 1 2J 



(9) 
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We use the following notations: 

k± = \ (h - h) , p = k 1 + k 2 , k' ± = ^(k' 1 -k' 2 ) , p' = k[ + k' 2 , (10) 

p p p'p' 

p2 p/2 / ± 1 M v !± M v 

* ~ S ' ~~ S 1 dfjLU ~ 9fMV ! 9uU ~ 9uU ) ) 

s s 

so one can write kj[ = k v g^ and k'^ = k' u g'J^ . In the centre-of-mass system, the integration 
can be re-written as 

'' d U^xA)^j^, (id 

4mj? 

where k! is the momentum of one of the quarks. The wave function reads: 

= ti S ' L ' J) (ki) . (12) 

Here Q^ 1 '. L .' l ?)(k±) is the moment operator for fermion-antifermion system [14] defined as 
follows: 

Q&\Ak±) = W^M) > ( 13 ) 
Q { ^lY\^)-^x^l a {k 1 _), 
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where 



k 2 



(2J- 1)!! 
J! 



h 1 - h 1 - h 1 - h 1 - h 1 - 



(14) 



2J 



~[ {ptiw^k^ ...k^ J + 9^ W3 k^ 2 k^ 4 . . . k^ + . . .) + 



A:1 



I "X ( fl -L fl X I.X 1.-L l.± _j_ 

(2 J — 1)(2J — 3) ^ww^wV w • • • V T 



IJ,i...lxj,a\ ro -i-J j2 Pi—tH-iM+i—PJ \ n, - L /y^ia 

2 V a 1 X( J -^ (k ) 

The interaction block can be expanded in a series in a full set of the t-channel operators: 

V (s, s', (k ± k' ± )) =J2Vi (s, s', (k ± k' ± )) I ®0 I , (15) 
1 

61 = I, 7 M , i^, ^75, 75 • 

The Bethe-Salpeter equation (8) is written in the momentum representation, it was solved 
in [4] in the momentum representation as well. Equation (8) allows one to use as an inter- 
action the instantaneous approximation or to take into account retardation effects. In the 
instantaneous approximation one has 

V (s, s', (k ± k' ± )) — V(t ± ), t ± = (k 1± - k' 1L )^-k 2± + k'^ (16) 

The retardation effects are taken into account, when the momentum transferred squared t 
in the interaction block depends on the time components of quark momenta (for more detail 
see Section 2.5 of [5] and the discussion in [13, 20, 22, 23]): 

V (s, s', (k ± k' ± )) — V(t), t = (k 1 - k'M-h + k' 2 ), . (17) 

It turned out that for the fit of the cc states we need only two interaction blocks with the 
following t-dependence: 



- (w^r ~ w^) ■ (18) 
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It also occurred that the results of the fit depended weakly on whether the t- or independence 
was used in (18). Because of that, for the sake of simplicity, we present below the results 
obtained in the instantaneous approximation substituting t — > t± in (18). 

Traditionally, the interaction of heavy quarks in the instantaneous approximation is pre- 
sented in terms of the potential V(r). The form of the potential can be obtained with the 
Fourier transform of (18) in the centre-of-mass system. Thus we have 

t± = -{k-k'f = -q\ 

/fW) = /(03 e "^MO, (19) 

that gives 

/(°) (r, n) = e-» r , (r, fj.) = r e^ r . (20) 

The potential used in [4] had the form 

V(r) = a + br + ce' lir , (21) 
where the constant and linear (confinement) terms read: 

a -> a J c (0) (r, n constant -> 0) , 

br -> b /^(r, fjL linear ->0) . (22) 

The limits ^constant , ^linear mean that in the fitting procedure the parameters ^constant 
and [innear are chosen to be small enough, of the order of 1-10 MeV. It ws checked that the 
solution for the considered states (n < 6) was stable when changing ^constant and fiimear in 
this interval, 1-10 MeV . 

In [4], charmonium wave functions were fitted in the following form: 

( ) 2 9 ( ) 

i=l 

Recall that k 2 is the relative quark momentum squared, s = 4m 2 c + 4k 2 , and n is the radial 
quantum number; (3=1 GeV~ 2 for all cc states. 

Two solutions with two types of the t-channel exchanges were used: 

Solution I : I <g) I , 7 M <g> 7 M , 75 <g> 75 , 

Solution II : I (g) I , 7^ <g> 7^ . 

Table 1 demonstrates the obtained values of the potential parameters (a, 6, c, /i) in Solutions 
I and II. In all solutions we put m c = 1.25 GeV. 

The measured masses of the cc-states and the results of the fit for n — 1,2, 3, 4, 5, 6 are 
displayed in Table 2. 
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Table 1: Parameters of the potential (in GeV units) 



Type of interaction (0/ (g> Oj) 


Solution 


a 


b 


c 


li 


Scalar (I ® I) 


I 


-1.527 


0.170 


1.013 


0.201 




II 


-1.417 


0.158 


0.883 


0.201 


Vector (7^ <g> 7^) 


I 


-1.539 





2.133 


0.401 




II 


-1.540 





2.130 


0.401 


Pseudoscalar (75 <g> 75) 


I 


-3.000 








0.201 




II 














Figure 5 shows the radiative decay transitions, which were included into fitting procedure 
in [4]. One can see the calculated numbers for partial widths and experimental values with 
errors used in the fit: the 20%-accuracy was accepted for the transitions ip{2S) — > 7x c j(lP) 
and 30%-one for x c j(1P) — > ^ipilS) (note that slightly smaller errors were obtained in the 
overall fit of Ref. [15]). 

The wave function parameters cf 1 ^ determined in (23) are presented in Table 3 for So- 
lutions I and II. Correspondingly, in Figs. 6 and 7 we demonstrate the wave functions for 
ip(nS) and rj c , Xco, Xci, Xc2- Comparing the wave functions depicted in Figs. 6 and 7, one can 
clearly see that Solutions I and II differ unsignificantly. We have carried out our calculations 
with two variants of the interaction in order to make clear that the description of cc system 
does not require a variety of the t-channel exchanges and inclusion of all the versions given by 
Eq. (15) would only result in the absence of convergence in fitting procedure. In particular, 
the considered example of two sorts of interactions demonstrate that the cc system does not 
require instanton-induced forces, which were needed for treating the mass splitting of ir, tj, 7/ 
[16]. 



3 Determination of the cc component of the photon 
wave function 



The vertex function of the transition 7 — > cc is represented with the following formula: 

<w«> = ± cww«) + 1+exp|( 4(, -,„)» • (24) 

where G^ nS )(s) = ^(„s)(s)(s - M| (nS) ) and M^ (n5) and ^(«s)( s ) are S iven in Tables 2 
and 3; C n , /3 7 and s are the parameters to be determined. 
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Table 2: Measured masses and results of the fit for Solutions I and II (in GeV units). 
Bold-faced numbers stand for masses which have been used in fitting procedure [4]. 



n 


ip{nS) 


r) c (nS) 


Xco(nP) 


XdinP) 


Xc2{nP) 




data 


fit 


data 


fit 


data 


fit 


data 


fit 


data 


fit 


i 


3.096 


3.1022 


2.979 


2.9776 


3.415 


3.3933 


3.510 


3.4962 


3.556 


3.5676 


2 


o.uoo 


O.O ( O ( 




q RO/lfi 
O.0Z40 




O.o4o0 




q anno 




q q/Iqc; 
o.y4yo 


Q 
O 


4.040 


4.0565 




4.0225 




4.1812 




4.2514 




4.3046 


4 


4.415 


4.3960 




4.3678 




4.5569 




4.6709 




4.8250 


5 




4.8465 




4.8233 




5.1185 




5.1886 




5.4758 


6 




5.4448 




5.4242 




5.7510 




5.8404 




6.2197 


1 


3.096 


3.1023 


2.979 


2.9772 


3.415 


3.3958 


3.510 


3.4979 


3.556 


3.5687 


2 


3.686 


3.6721 


3.594 


3.6238 




3.8447 




3.8957 




3.9434 


3 


4.040 


4.0470 




4.0139 




4.1705 




4.2406 




4.2976 


4 


4.415 


4.3801 




4.3527 




4.5488 




4.6605 




4.8188 


5 




4.8322 




4.8090 




5.1123 




5.1820 




5.4669 


6 




5.4336 




5.4135 




5.7460 




5.8379 




6.1986 



Table 3: Constants c- (in GeV units) for the wave functions of Eq.(23) in Solutions I, II 



i 


J/V(iS) 


V(2S) 


V(3S) 


V(4S) 


Vc(lS) 


Vc(2S) 


Xco(lP) 


Xco(2P) 


Xc2(lP) 


Xc2(2P) 


1 


8.15 


-18.87 


43.70 


-68.75 


-7.064 


-17.72 


27.17 


131.60 


-45.15 


-144.38 


2 


3.01 


9.16 


-271.18 


682.41 


0.08171 


4.60 


-12.77 


-480.42 


59.15 


569.17 


3 


-48.02 


266.00 


465.83 


-2356.65 


20.748 


247.04 


-100.08 


539.03 


141.47 


-587.98 


4 


76.13 


-729.46 


-92.87 


3873.67 


-14.849 


-653.28 


181.19 


-23.55 


-432.83 


-315.11 


5 


-49.16 


871.28 


-504.13 


-3421.29 


-22.992 


770.56 


-120.56 


-442.03 


481.62 


1126.02 


6 


9.17 


-566.22 


587.76 


1688.66 


41.115 


-500.42 


23.01 


410.77 


-286.31 


-978.20 


7 


4.81 


208.18 


-286.06 


-452.00 


-25.682 


185.63 


11.05 


-172.89 


96.19 


413.89 


8 


-2.65 


-40.78 


66.18 


57.34 


7.413 


-37.01 


-6.04 


36.23 


-17.21 


-87.74 


9 


0.3712 


3.320 


-5.948 


-2.183 


-0.8335 


3.094 


0.8269 


-3.081 


1.274 


7.454 


1 


8.14 


-19.41 


45.87 


-69.8 


-7.061 


-18.16 


27.79 


137.91 


-46.32 


-151.59 


2 


2.86 


12.56 


-295.78 


705.5 


0.29980 


7.16 


-15.64 


-520.67 


64.72 


624.91 


3 


-47.50 


260.19 


559.16 


-2482.2 


19.754 


243.34 


-95.29 


633.58 


132.65 


-753.22 


4 


75.58 


-729.85 


-261.10 


4172.7 


-13.001 


-655.51 


177.81 


-131.97 


-428.51 


-59.52 


5 


-48.97 


882.00 


-338.96 


-3793.5 


-24.954 


781.02 


-120.21 


-377.56 


484.48 


895.03 


6 


9.20 


-578.49 


494.77 


1949.0 


42.396 


-511.22 


23.98 


393.18 


-290.87 


-851.60 


7 


4.76 


214.44 


-256.26 


-554.9 


-26.199 


190.89 


10.43 


-172.64 


98.47 


372.46 


8 


-2.63 


-42.31 


61.177 


78.7 


7.527 


-38.26 


-5.88 


37.11 


-17.73 


-80.30 


9 


0.371 


3.468 


-5.6097 


-4.01 


-0.84433 


3.212 


0.8118 


-3.211 


1.320 


6.891 
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3.1 Decay ip{nS) — > e + e 

Partial width for the decay ip{nS) — > e + e~ reads: 



V(ip(nS) -> e+e - ) 



7T« 



^(n5)-^ /8 2 4 2 

M 5 \ 772 + o M V(n5) 



^+o M «)%SHeV(0) , (25) 



where a = e 2 /An = 1/137, \i e is the electron mass, and M^ nS ) is the charmonium mass. 
The transition amplitude i^( n s)-» e +e- (0), being determined by the process of Fig. 4, see [1], 
is equal to: 



z / — r ds / 

F^ nS )^e+e-(0) = %\ N c J J^^^nS) (s)G 7 ^ c5 (s) J 1 - 



4m2 /8 2 4 

£ -mi + -s 

s V3 c 3 



(26) 



2 — J r/s 

The wave function in (26) is normalized as follows: 

ds 

The coeficients Sf 1 given in Table 3 construct the wave functions obeying this normalization 
constraint. 



1 = / is^wW 1 



Ami /8 2 4 



S V3 mc + 3 S) • 



(27) 



3.2 Decay r] c (nP) — > 77 

Partial width for the decay r^ c — > 77 reads: 

rfac - 77) = ^« 2 MJ |i^ 77 (0)| 2 . (28) 
The transition amplitude is determined by the processes of Figs. la,b; it is equal to [1, 4, 17]: 



8 / — f ds 
F Vc{nP) _ 7J (0) = -\JN c m c J — * Tfc(nP) (s)* 7 _ >cg (s) 

4m? 



In 



(29) 



where iV c = 3. Recall that \l/ 7 ^ C c(s) = G 7 _> C c(s)/s . 

Normalization condition for pseudoscalar charmonium wave functions is as follows: 



ds 



Ami 



3^2 Vc(nP) 



S . 



(30) 



Coefficients presented in Table 3 give us \E' r?c ( n p)(s) obeying (30). 
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3.3 Decay Xco(nP) -> 77 

Partial width of the decay Xco — > 77 is equal to 



r(x c o -> 77) 



it a 



(o)l 2 , 



(31) 



L XcO 



with the quark transition amplitude (Figs. la,b) determined as follows [17, 18]: 

ds _ , , _ / \ I r~, —^7 „ 9 , v 7 ^ + \/s - 4m 2 . 



4rra? \ 

Normalization condition for scalar charmonium wave function reads: 



1 = 



4m? 



—^l , As)Jl-^(s-4m 2 c )m 2 c 

g n 2 XcO(nP)V I W s V c7 c 



^/i - ^/s - 4m 2 t 



(33) 



3.4 Decay Xc2(nP) -> 77 

Partial width x C 2 - *■ 77 is defined by two transition amplitudes: 

4 7Ttt 2 /l 



r(x C 2 -> 77) 



5M, 



Xc'2 



(34) 



which are determined by the processes of Figs. la,b [17, 19] and for the P-wave quark- 
antiquark states they read: equal to: 



7 {H) 
Xc2(nP)^77 



(0) 
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4m? 



Here 



J(°)( s ) = -2^5 (s - 4m 2 ) (l2m 2 + s) + 4m 2 (4m 2 + 3s) In 

j(2)( s ) = 4 V /S ^ ~ ^4 ( 5m 2 + 4m 2 ^ m 2 + s ) ln . 



s + y s (s — 4m 2 ) 
s — ^ ' s (s — 4m 2 ) 



s + Js (s — 4m 2 ) 



3 ' ' s — yj s (s — 4m 2 ) 

Normalization condition for the P-wave tensor cc system it is: 



00 

= 

7 16tt 2 

4m? 



dS T L(nP)(«)^Vl-^(8^ + 3 S )( S -4m 2 ) . 



(35) 



(36) 



(37) 
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3.5 The results of the fit 



By fitting the reactions involving the transition 7 — > cc, we have determined the parameters 
C n ,/3 7 , s defined in (24). For Solutions I and II, they are as follows (in GeV units): 





Solution I 


Solution II 


Ci 


= -4.945 


d = -4.995 


c 2 


= -2.893 


C 2 = -2.897 


c 3 


= -2.191 


C 3 = -2.179 


c 4 


= -2.300 


C 4 = -2.260 


c 5 


= -4.264 


C 5 = -4.368 


c e 


= -0.690 


C 6 = -0.479 




6 7 = 0.14 


6 7 = 0.15 




s = 11.9 


s = 22.5 



(38) 



Experimental values of partial widths included into fitting procedure versus those obtained 
in the fitting procedure are shown in Table 4. There are also predictions made for the 
two-photon decays of the first radial excitation states: ?? c (3594), xo(3847), %2(3947). 

Let us note that the decay Xc2(3556) — > 77 was not included into the fit because of 
controversity of the data. In the reaction pp — > 77, the value T(x2(3556) — > 77) = 0.32 ± 
0.080±0.055 keV was obtained [27], while in direct measurements, such as e + e~ annihilation, 
the width is much larger: 1.02 ± 0.40 ± 0.17 keV [24] , 1.76 ± 0.47 ± 0.40 keV [25] , 1.08 ± 
0.30 ±0.26 keV [26] . The compilation [21] provides the value close to that of [27]. The value 
found in our fit agrees with data reported by [24, 25, 26] and contradicts [27]. 

4 Conclusion 

We have carried out the calculations of radiative transitions, where the cc system participates, 
and compared the results with the experiment. The results are given in Table 5. In general, 
there is a good description of the data. Still, one should point to a disagreement for the 
following two cases: tp(2S) — > Xci0-P)l an d J/ip — > rj c (lS)^f. 

The calculation of partial width ip(2S) — > Xci0-P)l provides us with a value twice as 
large as given in [28, 29]. Such a disagreement may be related to either presumably much 
higher experimental error [28, 29] or a specific behaviour of the wave function of x cl (lP), 
that was not accounted for in [4]. 

Another discrepancy has been observed for the width of the transition J/ij} — > r] c {lS)~f. 
This is an Ml-transition, it is defined by the magnetic moment of the c-quark. One possibility 
to reduce the calculated value of J/ip — > r] c (lS)^ consists in the increase of the c-quark mass, 
another one is to include into consideration anomalous magnetic moment of c-quark. The 
hypothesis of the existence of anomalous magnetic moment at light quarks was suggested 
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Table 4: Calculated partial widths for nS and nP states versus experimental data (bold 
mass numbers stand for the predicted states) 



Decay 


T (keV) 


T exp (keV) 


IS: 


J/^(3096) - 


-> e + e~ 


5.444 (I) 
5.598 (II) 


5.40 ± 0.22 


2S : 


^(3686) -> 


e + e~ 


2.151 (I) 
2.210 (II) 


2.14 ± 0.21 


3S : 


^(4040) -> 


e + e~ 


0.756 (I) 
0.778 (II) 


0.75 ± 0.15 


4S : 


^(4415) -> 


e + e~ 


0.462 (I) 
0.498 (II) 


0.47 ± 0.10 


IS : 


^ c (2979) - 


_> 77 


6.979 (I) 
6.946 (II) 


7.0 ± 1.0 


2S : 


r? c (3594) - 


-> 77 


1.968 (I) 
1.034 (II) 




IP : 


Xco(3415) 


_> 77 


2.572 (I) 
2.440 (II) 


2.6 ± 0.5 


2P : 


Xco(3849) 


_» 77 


1.159 (I) 




3P : 


Xco(3845) 


_> 77 


1.021 (II) 




IP : 


Xc2 (3556) 


_> 77 


1.195 (I) 
1.155 (II) 




2P : 
3P : 


Xc2 (3950) 
Xc 2 (3943) 


_> 77 
_> 77 


2.051 (I) 
1.934 (II) 
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Table 5: Comparison of experimental data (in keV units) with our results and calculations 
of other groups. 



Decay 


Data 

X-/ Cat U ( A 1 


L CO LL1 LiQ 


LSCFlflOl 


LSfCMlOl 


RMCSIflll 


R MY Will 

±. U.VJ. \ V Mill 


NR ft 21 






26± 4 


22 


31-47 


26-31 


31 


32 


19.4 




XcH- 1 --^ J7 


25± 4 


59 


58-49 


63-50 


36 


48 


34.8 






20± 4 


19 


48-47 


51-49 


60 


35 


29.3 


V(2S) - 


»7c(15)7 


0.8±0.2 


0.4 


11-10 


10-13 


6 


1.3 


4.47 


Xco(lP) 




165±50 


169 


130-96 


143-110 


140 


119 


147 


Xcl(lP) 


- J/V(15)7 


295±90 


389 


390-399 


426-434 


250 


230 


287 




° I t r\- LkJ J I 


390±120 


229 


218-195 


240-218 


270 


347 


393 


J/V(15) 


- »?c(15) 7 


1.1±0.3 


4.1 


1.7-1.3 


1.7-1.4 


3.35 


2.66 


1.21 


J/>(1S) 


— » e + e~ 


5.4 ±0.22 


5.44 


5.26 


5.26 


8.05 


9.21 


12.2 


V(2S) -> 


e + e _ 


2.12±0.12 


2.15 


2.8-2.5 


2.9-2.7 


4.30 


5.87 


4.63 


V(3S) - 


e+e- 


0.75±0.15 


0.76 


2.0-1.6 


2.1-1.8 


3.05 


4.81 


3.20 


^(4S) -> 


e + e~ 


0.47±0.10 


0.46 


1.4-1.0 


1.6-1.3 


2.16 


3.95 


2.41 


»7c(15) - 


> 77 


7.0 ±0.9 


6.98 


6.2-6.3 


6.2-6.5 


4.2 


3.8 


19.1 


Xco(l^) 


_> 77 


2.6 ±0.5 


2.57 


1.6-1.8 


1.5-1.6 








Xc2(lP) 


77 


1.02±0.40±0.17[24] 
1.76±0.47±0.40[25] 
1.08±0.30±0.26[26] 
0.33±0.08±0.06[27] 


1.17 


0.3-0.4 


0.3-0.4 









rather long ago in connection with the description of the decay uj — > 7r°7 [30], see also 
discussion in [31]. 

In Table 5, the values of partial width are presented obtained by the other authors. 

In [10], the ideology of treating the cc system is similar to ours: the charmonium masses 
were fitted as well as the widths of radiative transitions. The results obtained in [10] depend 
on a chosen gauge for gluon exchange interaction — we demonstrate the results obtained 
both for Feynman (F) and Coulomb (C) gauges; different approaches used in [10] are reflected 
in the allowed accuracy intervals given in Table 5. 

In [11] the cc system was studied in terms of scalar (S) and vector (V) confinement forces 
- both variants are presented in Table 5. For the comparison we give in Table 5 the results 
obtained in the nonrelativistic approach to the cc system. 

Both in relativistic [10, 11] and nonrelativistic [32] approaches there is rather large dis- 
crepancy between the data and calculated values oiip(nS) — > e + e _ (in [10] the width of the 
transition J/ip — > e + e~ was not calculated but fixed). In our opinion, the reason is that in 
all above-mentioned papers, soft interaction of quarks was not accounted for — we mean the 
processes shown in Fig. 3b, c. In fact the necessity of taking into consideration the low-energy 
quark interaction, that is, the vector meson dominance in the transitions qq — > V — > 7, was 
understood decades ago but till now this procedure has not become commonly accepted even 
for light quarks: see, for example, [33, 34]. 
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Figure 1: a,b) Diagrams for the two-photon decay of cc state, c) Cuttings in the spectral 
integral representation. 




c 



(a) (b) 
Figure 2: Initial (a) and final (b) state interactions of quarks in the decay diagrams. 



(a) (b) (c) 

Figure 3: a) Graphical representation of the Bethe-Salpeter equation for the cc vertex; c,d) 
interaction of quarks in the vertex cc — > 7 and its approximation by the sum of transitions 
cc — > T,ip(nS) — > 7 . 
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Figure 4: Quark transition diagram for the process e + e — > ip . 




Figure 5: Radiative decays of the charmonium systems, which were taken into account in 
the fit [4]. The calculated decay partial widths are shown in the keV units (upper numbers) 
together with experimental data and errors accepted in the fit (lower numbers). 
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Figure 6: The cc wave functions for ip{nS). Solid and dashed lines stand for Solutions I and 
II. 
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Figure 7: The cc wave functions for r) c (nS), Xco(nP) and Xc2(nP). Solid and dashed lines 
correspond to Solutions I and II. 
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Figure 8: Vertex 7 — > cc ; solid and dashed lines correspond to Solutions I and II. 
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